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ANALYSIS OF FREQUENCY RESPONSE CHARACTERISTICS 
OF PROPELLANT VAPORIZATION* 
by Marcus F. Heidmann and Paul R. Wieber 
Lewis Research Center 

SUMMARY 

The dynamic response of a droplet vaporization process excited by traveling trans- 
verse acoustic oscillations is derived by a linear analysis. Results of a previous non- 
linear numerical study used to formulate the analytical model are reviewed, and a trans- 
fer function representative of the dynamics of the vaporization process is derived. The 
analysis provides dimensionless parameters related to propellant physical properties 
that characterize the dynamic behavior of the vaporization process. Application is made 
to the vaporization of heptane, oxygen, fluorine, ammonia, and hydrazine. The dynamic 
response of these propellants attains a peak value at a particular frequency. A compar- 
ison is made with a burning-rate process described by a characteristic time and an 
interaction index giving similar behavior. 


INTRODUCTION 

Acoustic mode instability in a rocket engine combustor occurs when the combustion 
energy is released in a manner that reinforces the acoustic oscillations. Various com- 
bustion processes that can limit or control energy release (chemical kinetics, drop 
burning, vaporization, jet breakup, drop shattering, etc. ) have been suggested as 
responsible for acoustic reinforcement. Many of these processes have been studied. 

For example, the dynamic behavior of several individual processes is reported in 
references 1 to 3, and process behavior measured in terms of combustor stability is 
given in references 4 to 6. Such studies have isolated specific problem areas in rocket 
engine instability that require additional study. 

Propellant vaporization is a process of particular interest. Reference 4 shows that 

^Presented at AIAA Second Propulsion Joint Specialist Conference, Colorado Springs, 
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variations in vaporization rate can significantly affect combustor stability. The vapori- 
zation process alone was studied more thoroughly in reference 7, in which the dynamic 
behavior of the process was obtained from a nonlinear numerical analysis by using a 
model that predicts the vaporization rate during the entire drop lifetime . Under certain 
conditions, the process exhibited dynamic behavior that may cause combustor instability. 
Such dynamic behavior was related to boundary conditions controlling the vaporization 
process of heptane drops. 

The nonlinear numerical analysis of reference 4 motivated the development of a 
linear dynamic analysis based on a simplified model for vaporization that gives the same 
dynamic behavior observed in reference 7. The advantage of the linear analysis is that 
it provides dimensionless parameters related to propellant properties that may be used 
to characterize and examine the dynamic behavior of the vaporization process for any 
propellant. 

In this report the results of the nonlinear numerical analysis of reference 7 are 
reviewed, and the linear analysis is presented and discussed. 


SYMBOLS 


2 

nozzle throat area, in. 

a^. acoustic velocity at nozzle throat, in. /sec 

b vapor pressure - liquid temperature coefficient, PjV^L 

C^Cg, constants 


Cg, etc. 

c specific heat of liquid, Btu/(lb)(°R) 

r 

p 

D molecular diffusion coefficient, in. 

f frequency, cps 

g gravitational constant, 32. 2 (lb mass)(ft)/(lb force)(sec^) 

k thermal conductivity, Btu/(in. )(sec)(°R) 

M mass of propellant being vaporized, lb 

J( molecular weight, (lb mass)/ (lb) (mole) 

N response factor 

Nu^ Nusselt number for heat transfer 

Nu, Nusselt number for mass transfer 

1X1 


2 


n interaction index 

o 

P pressure, lb/in. 

2 

P c combustion chamber pressure, lb/in. 

9 

Pj^ vapor pressure at propellant surface, lb/in. 

Pr Prandtl number 

q heat -transfer rate, arbitrary units 

q^ n heat transfer to propellant surface, Btu/sec 
q QU t heat transfer from propellant surface, Btu/sec 
R universal gas constant, 18 510 (in. -lb)/(°F)(lb)(mole) 
q initial drop radius, p 
Re Reynolds number 

r instantaneous drop radius, in. 

Sc Schmidt number 

s Laplace transform, d/dt 

T temperature of vapor film, °R 

combustion gas temperature, °R 
T l temperature of propellant, °R 

Tj. static temperature at nozzle throat, °R 

t time, sec 

t^, period of oscillation, 1/f, sec 

tgQ time to vaporize 50 percent of drop mass, sec 

Up final combustion gas velocity, ft/sec 

V volume, arbitrary units 

w vaporization rate, lb/sec 

Wj. nozzle mass flow rate, lb/sec 

Z correction factor for heat transfer 

a correction factor for mass transfer 

(8 vapor pressure - combustion chamber pressure parameter, dimensionless 

y ratio of specific heats 



Q 

p ^ gas density at nozzle throat, lb/in. 

X latent heat of vaporization, Btu/lb 

r characteristic time, sec 

r y mean drop lifetime, M/w, sec 

9 phase shift, deg 

ci> frequency, rad/sec 

Superscripts: 

’ denotes perturbation quantities (i.e. , x’ = (x - x)/x) 

— denotes mean values 


RESPONSE FACTOR 

A response factor can be defined that is one measure of the magnitude by which the 
combustion process can reinforce an acoustic oscillation. Such a response factor was 
introduced in the numerical analysis of reference 7 and will be used in this study to 
evaluate the dynamic behavior of the vaporization process. This response factor is 
based on the Rayleigh criterion for acoustic amplification by heat or mass addition. The 
Rayleigh criterion states that reinforcemnnt or amplification occurs when an excess of 
heat or mass is added while the pressure is greater than the mean value. A response 
factor taking into consideration the Rayleigh criterion is expressed for perturbations 
about a mean value by 


r v r t 

J j q'(V,t)P'(V,t)dtdV 

N = (1) 

yy [p’(v,t)j 2 dt dv 

where q’, the fractional heat or mass perturbation, is 

q* = l^ 
q 

and P 1 , the fractional pressure perturbation, is 
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p’ = 


p - p 


p 


When both q' 
finite volume, the 


and P' oscillate with the same periodicity and are uniform over a 
response factor is usually expressed for one period of oscillation as 


4 


f 



( 2 ) 


For sinusoidal oscillations in pressure, any heat release or mass flow process 
linearly related to pressure gives the following value of the integral: 


N = 


^max 

^max 


cos 6 


( 3 ) 


where 


P 'W = P max sin wt 


and 


q’(t) = q^ax sin ( wt + 

*> 

j In the analyses of propellant vaporization, the heat release rate q' is generally 

assumed synonymous with the mass release rate w*. This synonymity was assumed in 
the nonlinear analysis of reference 7. Also, an average value for the fractional perturba- 
tion in mass release rate over a finite volume was determined by a numerical integration 
(ref. 7). For these evaluations, w' varies nonlinearly with sinusoidal oscillations in 
pressure. The response factor as defined by equation (2), however, was approximated 
by numerical techniques over one period of the pressure oscillation. 
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Numerical values of the response factor within the range -1 to 1 were obtained in 
reference 7. Negative values indicate that excess heat is added when the pressure is 
less than the mean pressure (potential damping of acoustic oscillations), and positive 
values indicate that excess heat is added at pressures above the mean pressure (potential 
driving of acoustic oscillations). 

Some significance can be placed on actual numerical values of the response factor if 
the acoustic system is assumed to consist only of a heat or mass addition from vaporiza- 
tion (a potential acoustic gain) and a mass loss through an exhaust nozzle (a potential 
acoustic loss). If quasi -steady behavior is assumed, the combustor pressure and flow 
perturbations in a critical flow nozzle are in phase. The magnitude of these perturba- 
tions for adiabatic flow can be derived from 

w t = Vt a t 


where 



which gives 



If mass flow perturbations in the nozzle are considered analogous to heat or mass flow 
perturbations in the chamber, then, for sinusoidal oscillations in pressure, the response 
factor as given by equation (3) is 

N = _ fe 1 ) = (_0 ' 912, ’'= 1 - 2 

In a simple feedback loop corresponding to the assumed two-process system, the 
sum of such a negative nozzle response factor and the response factor of the vaporization 
process indicates whether an excess of mass is added when the pressure is higher or 
lower than the mean pressure, and thus, whether acoustic oscillations will decay or 
grow according to the Rayleigh criterion. 

The numerical result from this analysis of nozzle flow process suggests that the 
response factor for the vaporization process must exceed 0. 912 to cause combustion 
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instability. The exact dynamic behavior of nozzle flow and the magnitude of other 
acoustic losses and gains in rocket combustors is debatable, and therefore a precise 
criterion for stability cannot be established. In this study, therefore, attention will be 
confined to how the vaporization response factor varies with propellant properties and 
combustor operating conditions. 


NUMERICAL ANALYSIS 

The response factor as a function of oscillatory frequency was evaluated in the 
numerical analysis of reference 7 for a representative vaporization process and a 
specific acoustic mode of resonance. Calculations were made for n-heptane drops 
vaporizing in a cylindrical combustor containing heptane -oxygen combustion products. 
Pressure, velocity, and temperature oscillations associated with the first traveling 
transverse acoustic mode were superimposed on the normal combustion flow process. 
Drops of constant size were repetitively injected from positions uniformly distributed 
across the injector face. The drop vaporization theory developed in reference 8 was 
used. In these calculations, the acoustic oscillations affected the heat and mass transfer 
processes and the drop acceleration, which gave three-dimensional velocity components 
to the drops and caused perturbations in drop temperature and vaporization rate. 

The perturbation in vaporization rate at a given angular and radial position at which 
the pressure was uniform at an instant of time along the axis of the chamber was obtained 
from a summation of the vaporization rates of individual drops. A typical plot of the 
vaporization rate perturbation for times covering one complete pressure oscillation is 
shown in figure 1. Equation (4) was used to evaluate the response factor for each 
angular and radial position. The response factor for the entire chamber was taken as the 

average value of the response factor for all angular 
and radial positions. 

These calculations were made for a range of 
oscillatory frequencies and for a variety of boundary 
conditions affecting drop vaporization, including 
variations in combustor pressure, final combustion 
gas velocity, drop radius, initial drop temperature 
and velocity, and the amplitude of the pressure 
perturbation. The frequency response curves calcu- 
lated for these boundary conditions were correlated 
by a frequency parameter as shown in figure 2. An 
equally good single curve representation of results 
was obtained from a correlation based on drop life- 



Figure 1. - Nonlinear response of vaporiza- 
tion process (ref. 7). 
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Figure 2. - Correlation of response factor with vaporization parameters for heptane. 
Response factors are based on average rather than maximum pressure amplitude 
(ref. 7). 


time; this curve is given in figure 3, which shows the response factor as a function of a 
dimensionless time equal to the ratio of the drop half -lifetime (time to vaporize one -half 
the droplet mass) to the period of the oscillation. The response factors shown in 
figures 2 and 3 are larger than those reported in reference 7. The values have been 
recalculated and represent the average values of the response factors based on the 
pressure in each volume rather than the wall pressure used in reference 7. 

Figures 2 and 3 show the characteristic response of the vaporization process for 
heptane. For fixed boundary conditions, the response factor is a maximum at a specific 
frequency (corresponding to a dimensionless time of about 0.4). The response factor 
approaches zero at a lower frequency and a constant negative value at higher frequencies. 
The transition to negative response occurs at a dimensionless time of about 1. 4. From 
these response characteristics, the vaporization conditions conducive to driving and 
damping of acoustic oscillations can be readily established. It should be noted, however, 
that the maximum response factor (about 0. 8) is less than the nozzle loss (about 0. 912) 
previously calculated for simple dynamic behavior of exhaust flow. 

The factors responsible for this behavior of the vaporization process are of funda- 
mental interest. A description of the contributing factors was given in reference 7. 
Parameter groupings that can characterize specific behavior, however, are not readily 
deduced from a numerical analysis. Such groupings could be readily evaluated if a 
linear dynamic analysis were applicable to the problem being studied. 
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Figure 3. - Correlation of response factor with dimensionless time for heptane. Response 
factors based on average rather than maximum pressure amplitude (ref. 7). 


The primary restriction on a linear analysis is the dependency of the vaporization 
process on the absolute velocity difference between the drop and the combustion gases. 
Acoustic particle velocity is a maximum at both high and low pressures during the oscil- 
lation. For this reason, velocity difference and vaporization rate attain maximum values 
at conditions near both maximum and minimum pressure and cause the nonlinear behavior 
of vaporization rate shown in figure 1. An analysis of the perturbation curves in vapori- 
zation rate showed that the velocity difference contributions to the vaporization rate at 
high and low pressure were nearly equal and thus cancelled effects with regard to 
response factor evaluations. As an approximation, therefore, the velocity difference 
effects on vaporization rate can be ignored. With this assumption, a linearization of the 
vaporization equations is possible. The following section presents a linear dynamic 
analysis of a simplified vaporization process that is insensitive to perturbations in gas 
velocity. 
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LINEAR ANALYSIS 


A transfer function that characterizes the dynamic behavior of the vaporization 
process of various propellants will be derived by a linear analysis of a simplified vapori- 
zation model. 


Vaporization Model 

Analytical relations for the vaporizing drop presented in reference 8 will be used to 
develop linear equations for perturbations in the vaporization process of a single drop. 
The equations are assumed to apply to a drop in an array of drops of decreasing size 
down a combustion chamber. 


T 

V 


Drop mass: 


dM 

dt 


-w 


(5a) 


T « -W' 

V dt 


(5b) 


where 


r = (mean drop lifetime) 

w 


Vaporization rate: 


_ 2nDJ(r 
RT 


NuP, a 
m L 


where 


Nu m = 2 + 0.6(Sc) 1/3 (Re) 1/2 


and 



(6a) 

r 


(6b) 


(6c) 
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For this analysis, Nu is assumed proportional to (rPl 1 / 2 , 

m n c 

, and M is proportional to r , so that 


D 


is proportional to 


w = CjM^pJ/ 2 



(6d) 


and 


w' 


= I M’ + /3P* 
2 L 



(6e) 


where 


/3 = 



Vapor pressure: 


ln P L = C 2 ' 



(7a) 


P 


L 


= bT 


L 


(7b) 


where 


b = 


C 3 T L 


(T t . - C,)‘ 


Drop temperature: 


dT 


L 


dt 


1 

c pM 


^in %u0 


(8a) 
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If the variations in specific heat 
assumed negligible, then 


with drop temperature perturbations are 


CpT 

A 





’ q out 


(8b) 


where 


^in %ut 


= Aw 


Heat transfer to drop: 


qtn = 2tfkNu h r(T b - T L )Z 


(9a) 


where 


Nu h = 2 + 0. 6(Pr) 1 / 3 (Re) 1 / 2 


(9b) 


If the temperature difference (T, - T t ) and the correction factor for simultaneous 

D 1_J 1/2 

heat and mass transfer Z are constant, Nu^ is proportional to (rP c ) ' , and M is 
proportional to r 3 , then 


%, - C 5 MV2p J /2 


(9c) 


and 


qJ n = 1/2 M* + 1/2 P^ (9d) 

Heat transfer from drop: 

Vit = Xw ( 10a ) 

If variations in the heat of vaporization A with drop temperature perturbations are 
assumed negligible, then 


qout = 


w* 


(10 b) 


The model expressed by the linear equations may be described as follows. A drop 
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with a mass M and a mean lifetime r y is being vaporized. The mass of the drop is 
perturbating about an average value because the vaporization rate, which represents 
mass leaving the drop, is varying about an average level. The variations in vaporiza- 
tion rate depend on the mass of the drop and on the difference between the vapor pressure 
at the drop surface and the combustor pressure. The drop temperature determines the 
vapor pressure, with the temperature established in a heat reservoir. The heat to the 
reservoir depends on the combustor pressure and on the mass of the drop or the mass 
in the reservoir. Heat leaving the drop or reservoir is contained in the vaporizing 
propellant. 


Transfer Function 

The linear equations ((5b), (6e), (7b), (8b), (9d) and (10b)) can be combined by using 
the transform s = d/dt to give a transfer function involving the variables w' and p' : 


w 

Pi 


1 

2 


2t s 
v 


1 + 2r v s 


1 + (1 -2/3)^^r 
Aj3b ' 


i c p t l - 

1 + T v s 

A /3b v 


(ID 


This transfer function expresses the dynamic relation between w' and P . Solu- 

tions of the form w'(t) = w’ sin(cot + 0) can be obtained when P (t) = p' _ sin cot. 

max max 

The amplitude ratio and phase shift obtained from such solutions are 


r w^_ 

V P c/ 


2r v co 


max 


r 2I 1 / 2 

[l + (27 v co) 2 J 





1 + 

(1 - 2 P) CpTh T col 

2 


L A /3b v _ 


< 




\ X/3b V / 



1/2 


( 12 ) 


and 


c T c T 

8 = — - tan - '*' 2 t, co + tan - ^(l - 2/3) £ — r co - tan - * - j — t co 
2 v A/3b v A/3b v 


(13) 


The response factor defined by these solutions is given by the previously defined 
relation (eq. (3)) 
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TABLE I. - PROPELLANT PHYSICAL PROPERTIES AT EQUILIBRIUM 
DROP VAPORIZATION TEMPERATURE 


[Combustor pressure, 300 psi. ] 


Propellant 



Property 




Propellant 

Vapor 

Specific 

Latent 

Vapor 

Vapor 


temper- 

pressure at 

heat, 

heat of 

pressure - 

pressure - 


ature, 

propellant 

V 

vaporization, 

combustion 

liquid 


T L’ 

surface , 

Btu 


pressure 

temperature 


°R 

P L’ 

lb/in. 2 

(lb)(°R) 

Btu/lb 

parameter, 

P 

coefficient, 

b 

Heptane 

845 

133 

0.706 

93.8 

1.36 

8. 1 

Oxygen 

234 

275 

.421 

63.8 

4.43 

6. 5 

Fluorine 

220 

255 

. 376 

47.8 

2.99 

6.9 

Ammonia 

554 

205 

1.152 

483 

1.88 

8. 8 

Hydrazine 

859 

165 

.754 

1318 

1. 53 

10.0 



Dimensionless time, t 50 /t_ the dime nsionless time parameter . Response 

Figure 6. - Nonlinear dynamic response of factors for oxygen, fluorine, ammonia, and hydra - 
oxygen vaporization (ref. 7). 

zine, together with that for heptane, as functions of 
the dimensionless time parameter are given in figure 5. Physical properties used 
for the evaluations are given in table I. The listed properties are for the equilibrium 
drop temperature condition attained during steady vaporization at a combustor pressure 
of 300 psi (ref. 8). For comparison of the effects of physical properties on response, 
figure 6 shows the nonlinear evaluation of N for the condition of oxygen vaporization 
considered in reference 7. 
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DISCUSSION 


The response curve of heptane vaporization obtained by a linear analysis (fig. 4) has 
characteristic properties similar to those of the response curve obtained by the numeri- 
cal analysis (fig. 3). Differences in the response curves are primarily confined to the 
region of peak response. The linear analysis assumes that vaporization is initiated at an 
equilibrium drop temperature. A numerical evaluation in reference 7 for drops intro- 
duced at their equilibrium temperature gave a peak value of N = 0. 42, which is the 
value obtained by the linear analysis. The 800° R initial drop temperature calculation 
shown in figure 3 also gives comparable values. Heating of the drop from an initial 
injection temperature apparently introduces nonlinearity that increases the peak value of 
the response factor. 

A comparison of time bases for the response curves obtained by a linear analysis of 
a single drop and by a nonlinear analysis of a complete array of drops down the chamber 
shows that t^Q for this array is larger by a factor of about 4. 5 than for the single 
drop, where tgQ is the half -lifetime of the largest or initial drop in the array. The 
effective mean drop size of an array of drops is significantly less than the largest or 
initial drop. The relation between the largest drop and the effective mean drop depends 
on the distribution of drop sizes in the array. When equations (5a) and (6d) are used to 
specify time histories, however, an effective mean drop size of about one -fourth the 
initial drop size is obtained. Beyond the limitation of specifying an effective drop size, 
the linear analysis adequately describes the dimensionless times at which peak, zero, 
and negative response are obtained. These times are characterized by the parameter 

groupings appearing in the linearly derived transfer 
function (see eq. (11)). Figure 7 shows the contri- 
bution of these parameter groupings to the overall 
response of heptane vaporization. 

The term 2? v s/(l + 2T y s) characterizes the 
dynamics related to the quantity of propellant being 
vaporized. It establishes the response at low 
frequencies or small vaporization times. In this 
region, the quantity of propellant in the chamber 
varies inversely with the vaporization rate. In the 
extreme condition of zero dimensionless time, the 
propellant vaporizes as rapidly as it is introduced, 
and the response factor is zero. At frequencies 
equal to l/2r v and greater, the variations in the 
quantity of propellant being vaporized become small 
and do not restrict the vaporization rate. 



Figure 7. - Component characteristics of 
transfer function for heptane vaporization. 
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The parameter c T, t /XjSb characterizes the dynamics of heat storage within the 
P J-j V 

propellant. A frequency equal to l/(c p T jT^/X^b) essentially separates regions of oscil- 
lating and constant drop temperature. At lower frequencies the drop temperature oscil- 
lations approach an equilibrium condition with the oscillating environment. At higher 
frequencies, the temperature oscillations are reduced and eventually drop temperature 
and vapor pressure remain constant. 

A condition of constant vapor pressure gives the inverse effect of pressure oscilla- 
tions on vaporization rate from that obtained when vapor pressure is in equilibrium with 
the environment (eq. (6e)). With equilibrium conditions (generally low frequencies), an 
increase in combustor pressure increases the heat transfer to the drop, and the vapor 
pressure and the vaporization rate give a positive response factor. When vapor pres- 
sure is constant (generally high frequencies), an increase in combustor pressure 
suppresses vaporization rate, and a decrease in pressure accelerates the rate, a condi- 
tion comparable to flash vaporization. This process gives a negative response. The 
parameter (1 - 20)(cJT, r„/X/3b) characterizes the dynamics of these two interacting 
effects on vaporization rate. The values of j8 for the propellants evaluated in this study 
are greater than unity. Under such conditions, the transition from a vapor pressure to 
a combustion pressure controlled process occurs at a frequency that is lower by the 
factor l/(2/3 - 1) than the frequency characterizing the drop temperature behavior. 

These parameters, which characterize the dynamics of the vaporization process, 
are functions of propellant physical properties. From a knowledge of propellant physical 
properties, the dynamic response of various propellants can be estimated or computed. 
Oxygen vaporization (fig. 5) is one example. The region of positive response is broader, 
and the response factor attains a larger value than for heptane. The numerical evalua- 
tion (fig. 6), although incomplete, implies a similar behavior. 

Oxygen properties differ from those of heptane in several ways; however, the change 
in j3 (1. 36 for heptane; 4. 143 for oxygen) is the predominant factor causing a difference 
in response. The value of /3 for oxygen is large because the equilibrium vapor pressure 
for oxygen is a larger fraction of the combustor pressure than is the case for heptane. 
When /3 is large, the dimensionless time characterizing a constant vapor pressure in- 
creases. This phenomena increases and broadens the region of positive response. A 
large /3 also increases the perturbation in vaporization rate when vapor pressure is 
constant and thus increases the magnitudes of negative response as shown in figure 5. 

The comparison of various propellants (fig. 5) shows the region of positive response 
for hydrazine to extend over a broader region of dimensionless time than for heptane, 
or oxygen. Fluorine is nearly identical to oxygen. This comparison implies that a 
combustion process controlled by hydrazine vaporization would be unstable for a broad 
range of combustor designs. Reference 8 indicates, however, that decomposition of 
hydrazine could effect the vaporization process, which would also effect its dynamic 
behavior. 
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The dynamic response of other propellants for which physical properties are known 
or can be estimated can be evaluated in a similar manner. The linear analysis provides 
a convenient method of surveying the dynamic response of propellants and indicates the 
conditions in which the more precise numerical evaluation would be useful. 


STABILITY CRITERIA 

The dynamic response of the vaporization process does not, in itself, specify com- 
bustor stability. Stability criteria are acquired from a complete dynamic analysis of the 
overall combustion system. Such criteria have been derived in reference 9 for any 
combustion process that can be characterized by the flowing burning-rate expression: 

W = n[p^,(t) - Pj,(t - t)] (14) 


In this expression, the interaction index n and the characteristic time r specify 
a particular combustion process. The dynamic response for this burning rate expression 
is 


N = n(l - cos tu>) 


(15) 


A comparison for this dynamic response with the dynamic response of heptane vaporiza- 
tion is shown in figure 8. The region of positive response and particularly the peak 
response are of primary significance in establishing system stability limits. In this 
region, the dynamic response of the vaporization process can be approximated by the 
burning-rate expression (eq. (14)). For the comparison shown in figure 8, the vapori- 
zation process is characterized by n = 1/2 N max and r = 4. 5 r y . By a similar compar- 
ison with the correlation of the numerical results shown in figures 2 (p. 8) and 3 (p. 9), 

the value of r is approximated by 



T = t 5<)4 



3/2/ \l/3/ v 

/ 300 \ / 800 \ 

V p c/ V u f/ 


l/ 3 / xl 1 / 3 



(16) 


Figure 8. - Comparison of heptane vaporization 
with burning rate characterized by n and x. 


within the range of boundary conditions considered 
in the numerical analysis. The value of n is again 
equal to one -half the peak value of the response 
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factor, n ^ 0. 4. If the dynamic response of the vaporization process is characterized in 
this manner, the system stability criteria specified in terms of n and r are directly 
applicable to a vaporization controlled burning rate. For greater precision, a burning- 
rate expression may be postulated that more closely approximates the vaporization 
process and utilizes the method of analysis of reference 9 to establish stability criterion. 


CONCLUDING REMARKS 

The results of this study show that the dynamic response of the vaporization process 
can be approximated by a linear analysis that assumes the vaporization rate to be insensi- 
tive to gas velocity perturbations. The assumption is based on the net effect of velocity 
rather than actual insensitivity to velocity, that is, the changes in vaporization rate due 
to perturbation in velocity are the same at high and low pressures and, therefore, cancel 
with regard to dynamic response. The assumption appears valid for any transverse mode 
in the absence of any steady angular or radial flow within the cavity. With such steady 
flow, however, and in the case of longitudinal modes with axial flow of combustion gases, 
dynamic response is not insensitive to velocity perturbations (ref. 10). The linear 
analysis and the results of the numerical analysis do not necessarily apply to such condi- 
tions. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 1, 1966, 

128-31-06-02-22. 
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